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Abstract 



We show that the Hartree-Fock (HF) results cannot be reproduced within the 
framework of Kohn-Sham (KS) theory because the single-particle densities of finite 
systems obtained within the HF calculations are not v-representable, i.e., do not 
correspond to any ground state of a N non-interacting electron systems in a local 
external potential. For this reason, the KS theory, which finds a minimum on a 
different subset of all densities, can overestimate the ground state energy, as compared 
to the HF result. The discrepancy between the two approaches provides no grounds 
to assume that either the KS theory or the density functional theory suffers from 
internal contradictions. 
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The Hartree-Fock method (HF), first proposed in the pioneering works of Hartree and Fock [1,2] 
is known to be successful in calculating properties of electron systems, in particular, the ground state 
properties of atoms. Based on a variational principle, the HF method estimates the ground state 
energy E of an electron system from above, i.e., E HF > E, where E HF is the ground state energy 
calculated within the HF method. If the ground state wave function of N electrons is approximated 
by a single iV-electron Slater determinant, the HF solution delivers the minimum value Ehf on the 
set of all such determinants. Agreement, or otherwise, with the HF results is often used to estimate 
the success of other approximate computational schemes. 

The Density Functional Theory (DFT) exploits the one-to-one correspondence between the single- 
particle electron density and an external potential acting upon the system and relies on the existence 
of a universal functional F[p(r)] which can be minimized in order to find the ground state energy 
[3]. The Kohn-Sham (KS) theory goes further in reducing the problem of calculating ground state 
properties of a many-electron system in a local external single-particle potential to solving Hartree- 
like one-electron equations [3,4]. Successful solution of these equations allows to predict, at least in 
principle, the atomic, molecular, cluster and solid bodies binding energies, phonon spectra, activation 
barriers etc., see e.g. [5]. 

It is natural, therefore, to ask whether the HF ground state energy can be successfully reproduced 
in the Kohn-Sham approach. We note first that a universal density functional Fhf[p] can be defined 
with the help of the constrained-search technique [6]. Had the explicit form of F HF [p\ been available, 
the HF and the DFT approach would have yielded the same results for the ground state energy, 
E HF = E DHF , and the electron density, Phf(j) = Pdhf{ y )-> [6]. Here, E DHF and Pdhf{ v ) are the 
exact DFT Hartree-Fock energy and density respectively, while Phf( y ) is the HF density. Unfortu- 
nately, the correspondence theorem [3,6] establishes the existence of the functional only in principle, 
and provides no unique practical recipe for its construction. Rather, for practical calculations one has 
to resort to the KS approach. Exhaustive calculations [7-9] of the ground state energies of different 
atoms show that, if the KS approach is used, the resulting energy Ekshf usually exceeds the energy 
Ehf obtained by the HF method, 

Ershf > Ehf- (1) 

Recently, there have been suggestions that this disagreement may point to intrinsic flaws in both 
the DFT and the KS theories. One might suspect, therefore, that an exact local exchange potential 
does not exist for ground states of typical atoms (see [8,9] and references therein). We will, however, 
argue that the discrepancy (1) between the HF and KS is due to the different domains on which 
the respective functionals are defined. More specifically, we will show that while a KS density is 
^-representative, a HF density is not, i.e., it cannot be obtained as the ground state density of any N 
non- interacting electrons in a local potential. As a result, the KS method simply delivers a minimum 
on a different class of electron densities, and its disagreement with the HF approach does not indicate 
the existence of any internal contradictions either in KS or DFT approach. 

The purpose of this Letter is to show that HF results cannot be reproduced within the framework 
of KS theory because the single-particle densities of finite N electron systems obtained within the 
HF calculations are not -y-representable. This means that the HF densities cannot be obtained as the 
ground state density of any N non-interacting electrons in a local potential. Thus, the discrepancy 
between the two approaches, which manifests itself in the inequality (1), provides no grounds to 
assume that either the KS theory or the DFT suffers from internal contradictions. 
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We begin our study with considering the HF ground state energy which is given by the equation 
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We use an atomic system of units: e — m — h — 1, where e and m are electron charge and mass, 
respectively. Here N is the total number of electrons, rii are the occupation numbers: rii = 1 if the 
corresponding single-particle level is occupied, otherwise n» = 0. For an atom one has i>(r) = —Z/r, 
where Z is the nuclear charge. The exchange energy E x can be represented as follows 



1 r Xo(ri, r 2 , iiu) + 2vrp(r 1 )5(w)5(r 1 - r 2 ) rfw 



' x ~ 2 1 



|ri - r 2 | 
ki - r 2 | 



2tt 
rfridr 2 , 



dridr 2 



(3) 



and xo( r i ? r 2,^) is the linear response function, which is of the form 
Xo(ri, r 2 , W )=^ nj (l- n fc )^(ri)^ i (r 2 )0fc(r 2 )^fc(ri) ' ' 
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with u ik defined as uj ik = e k — £i, £k and functions 4>k{ T i) being respectively the one-particle energies 
and wave functions of the HF equations; and n is the infinitely small number, rj — > 0. Varying Eq. 
(2) with respect to the single-particle wave functions 0j(r), one obtains the HF system of equations: 
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These equations differ from an ordinary Schrodinger equation in two essential aspects: they are non- 
linear in 4>f F (r) and the second term under the sum on the left hand side of Eq. (5) that represents 
the so-called Fock's potential is non-local. 

The asymptotic behavior in r of 0^ F (r), contrary to the case of an ordinary one-particle 
Schrodinger equation, is not determined by Ef F and does not have the form 



J (r)| r _ >oo ~exp(- v /2|Ef ir |r). 



(6) 



On the contrary, it was shown in [10] that at r — > oo the function tpf F {r) is not determined by Ef F 
but behaves as 
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where Ep F is the energy of the so-called Fermi- level (with wave function <pp F (r)), which is the 
smallest binding energy of the occupied single-particle levels among all Ei in the considered system, 
and 
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Ci= £ / 0f F *(r)r'Y /m (r/r)0r(r)rfr, (8) 

m=—l 

and Y; m (r/r) is the I th order spherical polynomial. The uniform behavior of the occupied levels given 
by Eq. (7) is a consequence of the long range nature of the non-local Fock potential. The uniform 
behavior (7) leads to a number of very specific features of the ground state wave function, which, as 
we will see, cannot occur in the KS equations and makes it impossible for the KS equations (see Eq. 
(13)) to reproduce the HF density. That is, it makes the HF densities be non-f-representable. 

Consider next the HF calculations based on the constrained-search formulation of DFT which 
employs the HF density functional Fhf[p] [6]. The functional Fhf[p\ is obtained by minimizing the 
expectation value Fhf[p] = (^[p]\H\^f\p]) over all single-electron Slater determinants consistent with 
a fixed density p(r). The Hamiltonian H is of the form 

H = J- /^MV^rMr + \ / (9) 

From the manner in which the functional Fhf[p] is constructed it is clear that the resulting ground 
state energy E DHF equals that obtained in the HF method, E HF , 

E DHF = E H f = F HF [p DHF \ + J v{v)p DHF {v)dr. (10) 



Obviously, Pdhf(t) = Phf(t), where Phf(y) is the HF density obtained upon solving Eq. (5), while 
the determinant ty HF , which yields the minimum value, is composed of the eigenfunctions of Eq. (5) 
[6]. 

To give proof of the non-t> represent ability of the HF densities, we assume for a moment that 
Phf(?) is non-interacting t>-representable i.e., that it can be represented as the ground state density of 
N non-interacting electrons described by the Schrodinger equation with some local potential vl(t). 
Equivalently, we assume that the functional Fhf[p] is defined in the domain of f-representable 
densities. As a result, the HF functional Fjjf[p] can be represented as Fhf[p\ = F x [p] with F x [p] 
being a functional defined in the domain of f-representable densities, 

P[] T[] 1 v /- 0j(n)0?(r2)0i(ri)0.(r 2 ) f ( w . ( . 

F x [p\ = T k [p\ + - 2^ n 3 n i \ i i drirfr 2 + E x + v{r)p(r)dr. (11) 
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The second term in Eq. (11) being the Hartree term is obviously a functional of the density. The 
KS kinetic energy functional T k [p] is known to be defined in the domain of -y-representable densities 
[4,7] and so is the fourth terms in Eq. (11). Thus, the exchange energy E x given by Eq. (3) is also 
to be a functional E x [p\, defined on the -y-representable densities, as has been shown in [11], with the 
help of the representation (3). It has also been demonstrated in [11] that the variational derivative 
of E x [p] exists and can be evaluated explicitly to produce the KS exchange potential [11] 

(12) 

We can then proceed to obtain the eigenvalues an d the wave functions 0«(r) in Eqs. (4) and (11), 
by solving the KS single-particle equations 
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and compute the density p(r) as 

p(r) = 5>|<Mr)| 2 . (14) 

i 

We note that due to the constrained-search technique [6] the minimum of the functional F x [p] is 
given by 

F x [p] = (®ks\H\$ks) + / v(r)p(r)dr, (15) 

where &ks is a single iV-electron Slater determinant which delivers the lowest energy expectation 
value of H given by Eq. (9). This determinant is composed of the single-particle wave functions 0«(r) 
which are the solutions of one-particle equations (13). Thus, we are led to the conclusion that the 
determinant ^f HF has to coincide with the determinant ^ks- Then the HF wave functions 4>f F and 
the eigenvalues Ef F given by Eq. (5) must be equal to the wave-functions fa and the eigenvalues £j 
given by Eq. (13). However, it is seen from Eq. (7) that all single-particle HF functions have the 
same asymptotic behavior determined by the smallest orbital energy [10], even though the eigenvalues 
Ef F are not, in general, degenerate. On the other hand, eigenfunctions of Eq. (13) may exhibit such 
behavior only if the eigenvalues £j are degenerate as it follows from Eq. (6). As a result, we arrived at 
contradiction. This contradiction is resolved once we recognize that at least some of the HF densities 
are not non-interacting f-representable. Note that examples of non-f-representable densities were 
given in Ref. [13]. Obviously, a one-to-one correspondence between non-local potentials and local 
ones does not exist [14]. Therefore, if the explicit form the functionals F HF [p] and T k [p] was known, 
we would not have been able to obtain the HF ground state within the KS theory because they are 
defined on the different subsets of densities. For example, it is impossible to reproduce the HF single- 
particle eigenvalues Ef F within the KS theory [15]. In other words, the two approximate methods 
have different domains of applicability and are not amenable to a direct comparison. Consequently, 
the result (1) cannot be used to prove that the KS method or DFT is in any way deficient, as was 
suggested in Refs. [8,9]. 

A few remarks are in order here. A single A^-electron Slater determinant coinciding with &ks can 
be obtained in the optimized effective potential method with the local exchange potential Vopm [12]. 
Because of one-to-one correspondence that exists between the wave function, the density and the 
local single-particle potential [3,6], Vopm must coincide with V x given by Eq. (12), V x (r) = Vopm( v ) 
as has been demonstrated in [11]. Therefore, the ground state energy E OEP of a many-electron 
system calculated in OPM has to be equal to the corresponding energy E K shf calculated with the 
functional F x [p], Ekshf = Eoep- It is now tempting to assume that all three energies agree, i.e., that 
Ehf = Ekshf = Eoep [8,9]. But it is not the case since the HF densities are not f-representable. 
This fact is in agreement with numerical calculations showing that E HF < E KSHF [7-9]. One should 
also mention one special case where the HF and KS theories give the same answer. For a He atom, 
the HF potential acting on the occupied states is local and the HF density is f-representable, so 
that, Eoep = Ekshf = Eehf = Ehf- This observation is confirmed by numerical calculations [8]. 
For all other atoms, the HF potentials are non-local, and the HF densities are not f-representable. 
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Obviously, in that case, one has Eoep = Ekshf > Edhf = Ehf in accordance with the numerical 
calculations [8]. 

In summary, by clarifying the relationship between the non-local exchange HF potential and the 
local exchange KS potential, we have shown that the Hartree-Fock method cannot be reproduced 
within the framework of Kohn-Sham theory because the single-particle densities of finite systems 
obtained in Hartree-Fock calculations are not -y-representable. We have demonstrated that the fact 
that the KS calculations of finite electron systems lead to higher ground state energies cannot be 
used to infer the existence of inconsistencies in either KS or DFT theory. Most of the specific features 
of the HF method result from the non-local nature of the HF potential. For this reason, they provide 
no grounds to criticize the Kohn-Sham theory which deals with local single-particle potentials and 
v-representable densities. To conclude, it is worth mentioning that at present there is no compelling 
evidence to believe the HF method to be superior to the KS approach. Among the drawbacks of 
the HF theory is the well-known fact that the HF non-local single-particle potential acting on the 
unoccupied states falls off exponentially. As a result, the HF potential can only support very few 
unoccupied energy levels, which leads to difficulties in treating the excited states. By contrast, the 
KS theory does not suffer from the drawbacks inherent in the HF method. Thus, a further study is 
needed to clarify which type of behavior actually occurs in atoms with a large number of electrons 
and whether the failure to agree with the HF results can, indeed, be considered a fault of the KS 
theory. 
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